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Abstract 

Data from Type la supernovae, along with X-ray cluster estimates of the 
universal baryon fraction and Big Bang Nucleosynthesis (BBN) determina- 
tions of the baryon-to-photon ratio, are used to provide estimates of several 
global cosmological parameters at epochs near zero redshift. We show that 
our estimate of the present baryon density is in remarkably good agreement 
with that inferred from BBN at high redshift, provided the primordial abun- 
dance of deuterium is relatively low and the Universe is flat. We also compare 
these estimates to the baryon density at z ~ 1100 as inferred from the CMB 
angular power spectrum. 
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1 Introduction 



In the precision era of cosmology the accuracy in determining the key cos- 
mological parameters will be limited by our ability to constrain systematic 
errors and identify observations which break the many degeneracies between 
the global cosmological parameters and those related to specific models of 
inflation and/or structure formation. Redundancy can also provide valuable 
probes of unanticipated systematic uncertainties and may serve to separate 
global from model dependent parameters. In particular, although precision 
measurements of the anisotropy of the cosmic microwave background (CMB) 
promise statistically accurate determinations of many cosmological parame- 
ters, their interpretation may be limited by the extent to which parameter 
degeneracies can be broken and systematic uncertainties can be constrained. 
As a step in this direction, here we use current estimates of a restricted set 
of global cosmological parameters which are unaffected by "bias" (in mass 
versus light) and are independent of specific models of structure formation 
and specific theories of inflation to bound a variety of the key cosmological 
parameters. Such global cosmological constraints can then be employed in 
testing models of structure formation and inflation. 

The choice of which observations may provide the best constraints on the 
global cosmological parameters is time-dependent and subjective (see, for 
example, Steigman, Hata & Felten 1999 (SHF); Bahcall et al. 1999). It is 
our goal here to minimize the observational input while maximizing the pre- 
dicted output. To this end, we first utilize only the magnitude- redshift data 
from the type la supernovae (Perlmutter et al. 1997, 1999; Riess et al. 1997, 
1998; Schmidt et al. 1998; for a recent review and extensive references, see 
Filipenko & Riess 2000). In a FRW cosmology with an equation of state 
limited to two components, matter with zero pressure (p = 0) and vacuum 
energy (or a cosmological constant A) with negative pressure, p = —p, the 
assumption of flatness, in concert with the SNIa data, constrain the matter 
density (Qm) and the vacuum energy density (Qa) (Goobar & Perlmutter 
1995). Armed with fi M and we determine a variety of the other key cos- 
mological parameters such as the deceleration parameter (q = ^ — ^a) and 
the dimensionless age of the Universe (H t ). Furthermore, since a non-BBN, 
non-CMB constraint on the baryon density is of great current interest, we 
use estimates of the universal baryon fraction (/b) derived from X-ray obser- 
vations of galaxy clusters to bound the zero-redshift baryon density (White 
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et al. 1993 (WNEF); Steigman & Felten 1995; SHF). This independent de- 
termination of the baryon density (at zero redshift) may be compared to the 
high redshift estimates inferred from BBN (see Olive, Steigman, & Walker 
2000 and references therein) and from CMB anisotropy measurements (De 
Bernardis et al. 2000; Hanany et al. 2000; Lange et al. 2000; Jaffe et al. 2000 
and further references therein). Next, knowledge of Qm and Qb, along with 
the HST Key Project (Mould et al. 1999) constraint on the Hubble parame- 
ter permits us to bound other cosmological parameters such as the "shape" 
parameter T (Peacock & Dodds 1993; Sugiyama 1995; Peacock 1997; SHF) 
and to estimate the present age of the Universe to- 

We also explore a complementary approach by discarding the assumption 
of flatness (k — 0), and instead fixing the universal density of baryons from 
BBN. This along with estimates of the universal baryon fraction from the X- 
ray cluster data permits us to bound the total matter density, which may then 
be combined with the SNIa magnitude-redshift data to obtain constraints on 
the cosmological constant. With these constraints on Qm and Q\ we proceed 
to bound the 3-space curvature (f4 = 1 — (^m + ^a)), and the other global 
cosmological parameters such as go, H to (along with the age to), and T. 

2 Non-BBN, Non-CMB Cosmological Con- 
straints: Maximum Returns For Minimum 
Investment 

Two major groups have mounted systematic investigations of the high-redshift 
SNIa magnitude-redshift relation, the "Supernova Cosmology Project" (SCP) 
of Perlmutter et al. (1997, 1999) and the "High-Z Supernova Search Team" 
(HZT) of Schmidt et al. (1998). S. Jha and the HZT have kindly made avail- 
able to us the combined likelihoods and it is this joint data set we employ 
in our analysis. In our approach the computed quantities are the likelihood 
distributions for the various cosmological parameters (see Figures 2, 3, and 
6). Since none of our resulting distributions are perfectly gaussian, we report 
our quantitative results in two ways. We quote results in the form 
where A is the most likely value and the range from A — a_ to A + a + defines 
a 68% confidence region bounded by equiprobable points. We also quote the 
full 95% confidence range. 
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Figure 1: The 68% (solid) and 95% (dotted) contours in the Qa — plane 
allowed by the magnitude-redshift relation inferred from the joint (HZT and SCP) 
SNIa data. Closed universes lie in the region to the right and above the k = line 
while decelerating universes lie below and to the right of the qo = line (dashed). 

As is by now well known, the SNIa data identify a preferred region in the 
fi A — plane (see Figure 1) favoring an accelerating Universe (q < 0). The 
SNIa contours avoid Q\ = 0, and are cut by the flatness relation (k = 0): 
Qk = 1 — + = 0. If flatness is imposed as a constraint, the SNIa 
degeneracy between Qu and is broken and we find, in agreement with 
the SCP result of Perlmutter et al. (1999) and the HZT result of Filipenko 
& Riess (2000), 

Flat (k = 0) : Vl M = 1 - Vl A = 0.28±°$. (1) 

The likelihood distribution for is shown in panel a of Figure 2. The 
corresponding 95% range, 0.15 < Qm < 0.45, is entirely consistent with 
independent estimates (e.g., Cole et al. 1997; Carlberg et al. 1997; Bahcall 
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et al. 1999; Weinberg et al. 1999). 

Under the assumption of flatness the deceleration parameter is go = 



1 which, for the above value of f2 M , leads to q = — 0.58^q;}q. The likelihood 
distribution is shown in panel b of Figure 2; the corresponding 95% range is 
—0.35 > q > —0.77. The Universe is accelerating. 



1 

0.8 
\ 0.6 

i 

i 0.4 
0.2 





0.2 



0.4 0.6 

n.. 



0.8 



1 

0.8 
0.6 
0.4 
0.2 





-0.5 



0.5 





Figure 2: Likelihood distributions, normalized to unit maximum, for several global 
cosmological parameters as inferred from the SNIa data for a flat universe (k = 0). 
Panel a shows the total (baryonic plus CDM) density parameter J7m; panel b shows 
the deceleration parameter, qo; the age of the Universe in units of the Hubble age, 
Hgto, is shown in panel c; the shape parameter T is shown in panel d. 



A flat Universe has a dimensionless age (the age, to, in units of the "Hub- 
ble age", Ho" 1 ) of H t = ^sinh'Hg^) 172 - We find H t = 0.96t°;g 5 7 , 



in excellent agreement with the SCP result Hoto = 0.96lo;o7 (Perlmutter 
et al. 1998) and that from the HZT, H t = 0.94lg;^ (Riess et al. 1998; 
Filipenko & Riess 2000). The 95% range is 0.85 < H t < 1.13; the 
corresponding likelihood distribution is shown in panel c of Figure 2. For 



5 



H = 71 ± 6 kms _1 Mpc -1 , (HST Key Project: Mould et al. 1999) the age of 
the Universe is t = 13.2tJ;jj Gyr; the 95% range, 10.8 < t (Gyr) < 16.7, is 
in excellent agreement with that inferred from globular clusters (Chaboyer 
2000; Chaboyer & Krauss 2000). 

The shape parameter T (Peacock & Dodds 1993; Sugiyama 1995; Peacock 
1997; SHF) depends on Q M , tt B , h = Ho/lOOkms^Mpc- 1 , and the tilt 
parameter n. For k — and n = 1 we find T = 0.16^^1 and a 95% range 
of 0.08 < T < 0.26, which has considerable overlap with that inferred from 
observations of large scale structure (Fisher, Scharf, & Lahav 1994; Webster 
et al. 1998; Eisenstein & Zaldarriaga 1999; Efstathiou & Moody 2000). The 
likelihood distribution for T is shown in panel d of Figure 2. 

In summary, the SNIa data and the assumption of a flat Universe lead 
to a set of values for the global cosmological parameters Qm, ^a, and q , 
along with the large scale structure parameter T, which are consistent with 
other observational data and provide support for a Universe old enough to 
accomodate its oldest stars. 



2.1 X-Ray Cluster Baryon Fraction and the Baryon 
Density 

Rich clusters of galaxies, the largest collapsed systems in the Universe, pro- 
vide an ideal laboratory for exploring the universal baryon fraction (cf. 
WNEF; SHF). The X-rays observed from the hot intracluster gas may be 
used to estimate both the baryonic mass and the total mass of the cluster. 
Following WNEF, SHF have written the baryon fraction, / B , as 

/b = ^(1 + ^), (2) 

where /fie is the fraction of the total mass in the X-ray emitting hot gas, 
T is a baryon enhancement factor introduced by WNEF to account for the 
small offset between the universal and cluster baryon fractions, and the last 
term in eq. 2, taken from WNEF, accounts for the cluster baryons which 
are in stars rather than in the hot gas (while ignoring any possible contribu- 
tion from baryonic cluster dark matter). From a variety of hydrodynamical 
cosmological simulations Frenk et al. (1999) find for the offset between the 
universal and cluster baryon fractions: T = 0.92 ± 0.08. 
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SHF used Evrard's (1997) estimate of the hot gas fraction and chose 
f RG h 3/2 = 0.060 ±0.006. More recently, Mohr, Mathiesen and Evrard (1999) 
have summarized the observational results from an ensemble of clusters and 
also corrected for the effect of merger driven dumpiness (Mathiesen, Evrard 
& Mohr 1999) to find fnch 3 ^ 2 = 0.075. Although their quoted formal statis- 
tical uncertainty is small, systematic errors likely dominate the error budget. 
For example, while the latter estimate has been obtained using the "isother- 
mal beta model" as a total mass estimator, Evrard (Private Communication), 
preferring the virial theorem mass estimator, finds /hg^ 2 — 0.056. In an at- 
tempt to account for this uncertainty, here we adopt the average of the Evrard 
(1997) and the Mathiesen, Evrard, & Mohr (1999) determinations, along 
with a correspondingly generous error estimate: fnch 3 ^ 2 = 0.066 ±0.013. As 
gravitational lensing observations of X-ray clusters improve (see, for example, 
Tyson, Kochanski & Dell'Antonio 1998) this uncertainty should be reduced 
considerably. Observations of the Sunyaev-Zeldovich effect in clusters also 
promise more accurate determinations of the cluster hot gas fraction (see, for 
example, Grego et al. 2000). 

The present ratio of baryons to (CMB) photons is parameterized by r]i = 
10 10 n B /n y which, for a present CMB temperature of 2.725 K (Mather et 
al. 1999), may be written in terms of the baryon density parameter (£l B ) and 
the Hubble parameter (h) as, 

Vl0 = 274 n B h 2 = 274 (f B h 2 ) Q M . (3) 

Combining the above X-ray cluster estimates with the HST Key Project 
determined Hubble parameter (Mould et al. 1999), we obtain the distribution 
for f B h 2 . We find f B h 2 = 0.065±g:8ii and the 95% range is 0.037 < f B h 2 < 
0.099. Convolving this with the previously determined distribution (SNIa, 
k — 0) for Qm we find, 

Flat(k = 0): ^o = 4.8H;l (Q B h 2 = 0.0181°;°°?). (4) 

The likelihood distribution for rj is shown in Figure 3. While the 95% 
range, 2.1 < r] 10 < 9.1 (0.008 < Vl B h 2 < 0.033), is very broad, this de- 
termination of the baryon density (at zero redshift) is independent of, and 
complementary to, those from BBN (z ~ 10 8 ) and the CMB (z 10 3 ). As 
will become clear in the next section (see, also, Fig. 3), this estimate of the 
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Figure 3: The likelihood distributions, normalized to unit maximum, for the 
baryon-to-photon ratio r/io = 27AD,^h 2 . The solid curve is for the low-D BBN 
case; the dashed curve is for the high-D BBN case (see Sec. 3); the dotted curve is 
for the non-BBN case (SNIa h k = 0); the dot-dashed curve is the CMB inferred 
range for Q B h 2 = 0.032 ± 0.005 from Jaffe et al. (2000). 

baryon density is in excellent agreement with that determined from BBN us- 
ing the "low" deuterium abundance of Buries & Tytler (1998a,b; hereafter, 
BT). Although offset from the value of the baryon density consistent with 
BBN and the "high" deuterium abundance of Webb et al. (1997), given the 
large errors for both determinations there is still considerable overlap. 

Fukugita, Hogan, and Peebles (1998) have attempted to inventory baryons 
at z w 0, finding a range (for h = 0.70) 0.007 < n B < 0.041. Although this 
range (corresponding to 1.0 < 7710 ^ 5.5) has considerable overlap with ours, 
there is a hint that some dark baryons may have escaped the Fukugita, 
Hogan, & Peebles (1998) inventory. 

Measurements of the "large V multipoles of the CMB temperature angu- 
lar power spectrum can also probe the cosmological parameters. The recent 
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Boomerang (De Bernardis et al. 2000; Lange et al. 2000) and Maxima-1 
(Hanany et al. 2000) CMB anisotropy results point to a "low" second acous- 
tic peak. In the attempts to account for a second peak which is low relative 
to the first peak in the CMB anisotropy spectrum, there is some degener- 
acy between "global" cosmological parameters such as the baryon density, 
and model dependent parameters such as "tilt". Nevertheless, all multi- 
parameter fits to the Boomerang and Maxima data appear to point to a 
"high" baryon density. For example, Jaffe et al. (2000) find for their flat 
model (k = 0) which provides the best fit to the combined Boomerang, Max- 
ima and COBE DMR data, Vt B h 2 = 0.032 ± 0.005. On the assumption of 
gaussian errors, (almost certainly wrong) the likelihood distribution for this 
result is plotted in Figure 3; note the significant overlap between the baryon 
density determined at present from the SNIa data and that from the CMB. 
It is a remarkable confirmation of the (current) standard model of cosmology 
that these two, independent determinations of the baryon density, "mea- 
sured" at vastly different epochs in the evolution of the Universe, agree so 
well. Given the large uncertainties in each determination, the hint of a CMB 
challenge to the zero-redshift baryon density is only minor at present. 



3 Using The BBN-Inferred Baryon Density 

In this section we drop the assumption of flatness and replace it with es- 
timates of the baryon density from BBN. Again we use the X-ray cluster 
inferred baryon fraction for an estimate of the universal baryon fraction, but 
this time we combine /b with the BBN measure of f^B^ 2 (utilizing the esti- 
mates of the primordial abundance of deuterium) to obtain an estimate of 
the overall matter density (see Figures 4 & 5): 

= ^ = 0.057(1 ± 0.22)r7io, (5) 

JB 

where we have accounted for the uncertainties in /b and h as discussed above. 
Armed with the SNIa data and this BBN-related estimate of Qu, we proceed 
to determine the likelihood distributions for the global cosmological param- 
eters (Q M , ^a, fifc, go, r, H t ). 

It has long been known that deuterium provides an ideal baryometer 
(Reeves, Audouze, Fowler and Schramm 1976). Because deuterium has only 
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Figure 4: Same as Figure 1 with the 68% and 95% contours for Qm inferred from 
BBN and the X-ray cluster baryon fraction. The high-D, low-r/ constraint bands 
(at the left) are in red/dashed, while the low-D, high-?] bands (at the right) are in 
blue/dotted. 

been destroyed in the course of galactic evolution since the epoch of primor- 
dial nucleosynthesis (Epstein, Lattimer and Schramm 1976), any deuterium, 
observed anywhere in the Universe, provides a lower limit to its primordial 
abundance. In the high-redshift, low-metallicity environment of the QSO 
absorption-line systems the bulk of the gas is unlikely to have been processed 
through stars, so that the observed deuterium abundance is likely to be very 
nearly the primordial value. For two such systems BT derive a statistically 
accurate estimate for the primordial-D abundance: (D/H)p = 3.4±0.3 x 10~ 5 . 
The ~ 8% observational error is well matched to the comparable error in the 
BBN prediction (Hata et al. 1996). Due to the steep dependence of the BBN- 
predicted D-abundance on the baryon-to-photon ratio r], a 10% uncertainty 
in (D/H)p results in a baryon abundance determined to 6%. Combining the 
observational and theoretical errors in quadrature, the BT data lead to an es- 
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timate of the baryon abundance accurate to 7%. The likelihood distribution 
for T]io for this low-D case is shown in Figure 3 for 

Low-D: r7io = 5.1 ±0.4 (tt B h 2 = 0.019 ± 0.001). (6) 

Note the remarkable agreement between this BBN-determined baryon den- 
sity and that obtained in the previous section using the X-ray cluster baryon 
fraction and the SNIa and k = estimate (see Figure 3 for the correspond- 
ing likelihood distribution). The baryon abundance in the 10 - 15 Gyr old 
universe and that inferred for when the universe was only a few minutes old 
are virtually identical (as they should be!). 




Figure 5: Same as Figure 4, but showing the 68% and 95% contours for the joint 
fit of the SNIa data with the BBN/X-ray cluster determined values of Qy[. 

In contrast, a higher primordial deuterium abundance, (D/H)p = 20 ± 
10 x 10 -5 is inferred from HST observations of a different, lower redshift QSO 
absorption-line system (Webb et al. 1997). While the high deuterium inferred 
for this absorption system may well be contaminated by a low column density 
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hydrogen "interloper" , it could be that this higher abundance is primordial. 
If it were, then that found by BT would be low because the gas in their 
absorbing systems had been cycled through stars where the deuterium was 
destroyed. It would be surprising that the two high redshift systems are 
"evolved" while the one lower redshift system is not (or, is less evolved). 
Clearly more data will be crucial in resolving this question. Here we will 
attempt to use consistency with other cosmological parameters as a potential 
way to discriminate between the two possibilities ("low-D" and "high-D"). If 
the Webb et al. (1997) deuterium abundance is adopted, the corresponding 
baryon abundance is both lower and somewhat less accurately constrained 

High-D: r] W = 1.7 ±0.3 (Vt B h 2 = 0.006 ± 0.001). (7) 

From Figure 3 it is clear that the overlap with the non-BBN estimate from 
the previous section is poor. 

3.1 Global parameters from low-D 

Here we adopt the BBN-determined baryon density corresponding to low- 
D and we use the X-ray cluster baryon fraction to estimate the total mass 
density. The corresponding 68% (95%) bands are shown in Figure 4 and 
the corresponding contours in the — Qu plane are shown in Figure 5. 
Notice that the low-D determined contours are completely consistent with 
a flat Universe. The likelihood distributions for several global cosmological 
parameters are shown by the solid curves in panels a - / of Figure 6. 

For the total (baryonic plus CDM) mass density we find fi M = 0.26lg;og; 
the 95% range is 0.16 < Qu < 0.48. This is in excellent agreement with 
the k = identified range from the SNIa data and with other, independent 
determinations of the mass density (Cole et al. 1997; Weinberg et al. 1999) 
providing support for BBN with low-D. 

In this section we have dropped the restriction to a flat universe, thereby 
gaining an independent determination of the vacuum density (cosmological 
constant). The likelihood distribution for Q\ is shown by the solid curve 
in panel b of Figure 6. We find strong support for a non-zero value: Q\ = 
0.741^; the 95% range is 0.35 < A < 1.11. 

Having removed the constraint of flatness we can now ask if low-D BBN 
is consistent with a flat Universe. It is clear from the solid curve in panel 
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c of Figure 6 that a flat universe (Qk = 0) is completely consistent with 
the matter density determined from low-D BBN and the SNIa magnitude- 
redshift relation. Indeed, our best fit is for Q k = O.OOlg;^ and the 95% range 
is -0.53 < tt k < +0.44. 

Given the similarity between the non-BBN parameters and those identi- 
fied by low-D BBN, it is not suprising that here, too, we confirm an acceler- 
ating universe. The likelihood distribution for q is shown by the solid curve 
in panel d of Figure 6. Our best fit is for q = — 0.60l ;]J and the 95% range 
is -0.24 > g > -0.91. 
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Figure 6: Likelihood distributions for several global cosmological parameters de- 
termined from the joint fit of the SNIa data, the X-ray cluster baryon fraction, and 
the BBN constraints on the baryon density from low-D (solid curves) and high- 
D (dashed curves). Panel a shows the total mass density parameter Qm', panel 
b shows the density parameter associated with a cosmological constant Q\; the 
"flatness" parameter = 1 — (J7m + ^a) is shown in panel c; the deceleration 
parameter qo is shown in panel d; in panel e is the dimensionless age Hoto; the 
shape parameter T is shown in panel /. 
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The dimensionless age of the universe, Hoto, is shown in panel e of Figure 
6. For this low-D universe we find H to = 0.97io!o5! with a 95% range of 
0.88 < H t < 1.10. Adopting the HST Key project value of H = 71 ± 6 
kms _1 Mpc _1 (Mould et al. 1999), we find for the age of the Universe to = 
13.31^ Gyr with a 95% range extending from 11.0 Gyr to 16.5 Gyr. 

Finally, we may predict the value of the shape parameter corresponding 
to this low-D case. The likelihood distribution is shown in panel / of Figure 
6. It corresponds to a somewhat low value T = 0.15io!o4> but the 95% range, 
0.08 < T < 0.28 is perfectly consistent with observational estimates (Peacock 
k Dodds 1994; Fisher, Scharf & Lahav 1994; Maddox, Efstathiou & Suther- 
land 1996; Webster et al. 1998; Eisenstein & Zaldarriaga 1999; Efstathiou & 
Moody 2000). 

Very recently O'Meara et al. (2000) have reported the detection of deu- 
terium in another high redshift, low metallicity QSO absorption system. The 
deuterium abundance they derive, D/H = 2.54 ± 0.23 x 10~ 5 , is some 3a 
smaller than the previous mean D/H for the "low-D" absorbers. Accounting 
for the dispersion among the three, low-D absorbers, O'Meara et al. sug- 
gest a revised primordial deuterium abundance: D/H = 3.0 ± 0.4 x 10" 5 , 
which for BBN corresponds to a baryon-to-photon ratio of r] W = 5.6 ± 0.5, 
or fish 2 = 0.020 ±0.002. This small shift in the baryon density has minimal 
effect on the quantitative conclusions reached above. For example, while the 
95% ranges for Q M and Q\ shift upwards by approximately 0.03 and 0.04 
respectively, that for Q k moves down by « 0.05, and the ranges of q , T, 
H t , and t are virtually unchanged. 

In summary, the low-D BBN estimate of the baryon density, in combina- 
tion with the X-ray cluster estimate of the baryon fraction and the SNIa data 
favors a flat (k — 0), accelerating (g < 0) Universe whose age (to > 11 Gyr) 
is consistent with estimates of the ages of the oldest stars and with a shape 
parameter (0.1 < T < 0.3) in agreement with estimates from observations of 
large scale structure. 

3.2 Global parameters from high-D 

If, in contrast to the above, the BBN-determined baryon density correspond- 
ing to high-D (Webb et al. 2000) is used, we find a very low matter (baryons 
plus CDM) density, fi M = 0.09 ± 0.03; the 95% range is 0.04 < fi M < 0.17. 
The corresponding vertical band for Qm is shown in Figure 4, the related 



14 



— contours are shown in Figure 5, and the likelihood distribution for 
Qm is the dashed curve in panel a of Figure 6. Such a low value for the 
total mass density is in conflict with independent determinations (Bahcall 
et al. 1999) and with the comparison between models of structure forma- 
tion and observations of large scale structure (Cole et al. 1997; Weinberg et 
al. 1999). 

The corresponding likelihood distributions of the other global cosmolog- 
ical parameters are shown in Figure 6. It is clear from panel c that a flat 
Universe is very unlikely (95% range: 0.08 < Qk < 0.78); Qk < is excluded 
at the 99.3% confidence level if the high-D BBN-determined baryon density 
is assumed to provide a fair estimate of the present baryon density. 

It is clear from Figure 6 that while the high-D, BBN determined values of 
f^A, %, and H t differ little from their low-D counterparts (at 95%: 0.15 < 
tt A < 0.79; -0.10 > g > -0.74; 0.92 < H t < 1.19), the corresponding 
range for T is very low (0.01 < T < 0.10), in conflict with observational 
estimates (Peacock & Dodds 1994; Fisher, Scharf & Lahav 1994; Maddox, 
Efstathiou & Sutherland 1996; Webster et al. 1998; Eisenstein & Zaldarriaga 
1999; Efstathiou & Moody 2000). 

In summary, in contrast to the low-D case considered above, the high-D 
estimate of the baryon density in combination with the X-ray cluster estimate 
of the baryon fraction and the SNIa data does not favor a flat (k = 0) Universe 
and prefers a low matter density (Qm < 0.17) and a small value of the shape 
parameter (r < 0.10) which are in marginal conflict with estimates from 
observations of large scale structure. These latter data argue against high-D 
and favor the low-D estimate of the BBN-predicted primordial deuterium 
abundance. 

4 Discussion 

Data from observations of the SNIa magnitude-redshift relation identify a 
region in the Q\ — plane consistent with a flat universe (k — 0), which 
disfavors a vanishing cosmological constant (Q\ ^ 0), and points towards 
an accelerating expansion (qo < 0) of the present universe. Imposing the 
assumption of flatness (k = 0) breaks the degeneracy between and fi M , 
identifying preferred values for each parameter separately, leading to well- 
defined predictions for a variety of other global cosmological parameters (see 
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Figures 1 & 2). In particular, we have supplemented the SNIa data with 
estimates of the cluster baryon fraction to relate the total mass density to 
the baryonic mass density, finding a value (determined at redshifts < 1) 
consistent with that suggested by BBN (fixed at redshifts > 10 8 ) for the 
case of a low primordial abundance of deuterium (r] 10 pa 5.1, Vt^h 2 ~ 0.019). 
Although the CMB-inferred value of the baryon density ( Jaffe et al. 2000) has 
a large uncertainty, and the result is dependent on non-global cosmological 
parameters, it too is roughly consistent with this value (albeit on the high 
side by some 60%). In contrast, the BBN-determined value of the baryon 
density in the case of a high primordial abundance of deuterium (1710 ~ 1.7; 
Q B h 2 pa 0.006) is inconsistent with both the SNIa - X-ray cluster result and 
that from the CMB anisotropy measurements. Our results are summarized 
in Figure 3 where we show the likelihood distributions for for the cases 
considered here, as well as the estimate from the CMB anisotropy (Jaffe et 
al. 2000). It is clear that there is excellent overlap between the non-BBN 
SNIa range (assuming flatness) and the low-D BBN range. In contrast, high- 
D BBN combined with the X-ray cluster baryon fraction leads to a very low 
estimate of the total matter density which, when combined with the SNIa 
data, appears inconsistent with a flat universe. 

In all cases considered here (non-BBN, low-D BBN, high-D BBN) the 
SNIa data strongly favor an accelerating Universe (go < 0). For all cases 
the "age problem" has evaporated (H to ~ 1, to ~ 13 — 14 Gyr). However, 
although both the non-BBN and the low-D BBN results predict a value 
for the shape parameter (r pa 0.2) consistent with observations (Peacock & 
Dodds 1994; Fisher, Scharf & Lahav 1994; Maddox, Efstathiou & Sutherland 
1996; Webster et al. 1998), the high-D prediction (r « 0.05) is low. All in 
all, our analysis supports a flat Universe (Qk ~ 0), presently dominated by 
a non-zero cosmological constant (JIa pa 1 — il M pa 0.7) or some other form 
of "dark energy" , and a baryon density consistent with that inferred from a 
comparison of BBN and the low value of the deuterium abundance found by 
BT (7/10 ~ 5; VL-qK 1 pa 0.02). While this baryon density is roughly a factor of 
3 larger than that inferred from BBN and the high deuterium value of Webb 
et al. (1997), it is within less than a factor of two from that derived from 
the CMB anisotropy data (Jaffe et al. 2000). Indeed, agreement between the 
baryon abundance at epochs in the evolution of the Universe (BBN, CMB, 
z pa 0) separated by some 10 - 14 orders of magnitude in time provides 
powerful support for the standard hot big bang cosmological model. 
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